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We present measurements of the branching fractions and charge asymmetries of B decays to 
all D^*'D^*' modes. Using 232 million BB pairs recorded on the T(4S') resonance by the BABAR 
detector at the e+e~ asymmetric B factory PEP-II at the Stanford Linear Accelerator Center, we 
measure the branching fractions 



B{B° ^ D'+D'-) = (8.1±0.6±1.0) X 10"^ 
B{B" ^ 0*^0"^) = (5.7±0.7±0.7) X 10"^ 



B{B°- 


-^ D+D- 


B{B+ ^ 


D'+D"' 


B{B+ - 


-. D'+W 


B{B+ - 


-. D+D"' 


and B(B+ 


^ D+W 



= (2.8 ±0.4 ±0.5) X 10" 

= (8.1 ±1.2 ±1.2) X 10" 

= (3.6 ±0.5 ±0.4) X 10" 

= (6.3 ±1.4 ±1.0) X 10" 

= (3.8 ±0.6 ±0.5) X 10" 

where in each case the first uncertainty is statistical and the second systematic. We also determine 
the limits 

B{B° ^ D'^D*") < 0.9 X 10"", 

B{B° -^ D*"d") < 2.9 X 10'"^ 

and Z?(B" -^ D^D") < 0.6 x 10"", 

each at 90% confidence level. All decays above denote either member of a charge conjugate pair. 
We also determine the CP- violating charge asymmetries 



A{B° ^ D'^D"^) = 


= 0.03 ±0.10 ±0.02, 


A{B+ -> D*+^*°) = 


= -0.15 ±0.11 ±0.02 


A{B+ -^ D*+^") = 


= -0.06 ±0.13 ±0.02 


A{B+ -^ D+D"^) = 


= 0.13 ±0.18 ±0.04, 


and A{B+ -^ D+W) = 


= -0.13 ±0.14 ±0.02 



Additionally, when we combine these results with information from time-dependent CP asymmetries 
in B" — > !)(*'+ £)(*'" decays and world-averaged branching fractions of B decays to Ds D^*' modes, 
we find the CKM phase 7 is favored to lie in the range [0.07— 2.77] radians (with a ±0 or +n radians 
ambiguity) at 68% confidence level. 

PACS numbers: 13.25.Hw, 12.15.Hh, ll.30.Er 



I. INTRODUCTION 

We report on measurements of branching fractions of 
neutral and charged i?-meson decays to the ten double- 
charm final states D^*^D'-*\ For the four charged B de- 
cays to _d(*)d(*) and for neutral B decays to D*^D^, 
we also measure the direct CP-violating time-integrated 
charge asymmetry 



A 



T+ 



CP 



T- 



^^ 



(1) 



where in the case of the charged B decays, the superscript 
on r corresponds to the sign of the B^ meson, and for 
D*±Z)T, r+ refers to D*-D+ and T" to D*+D-. 

In the neutral B -^ £)(*)+£)(*)- decays, the in- 
terference of the dominant tree diagram (see Fig. la) 
with the neutral B mixing diagram is sensitive to 
the Cabibbo-Kobayashi-Maskawa (CKM) phase f3 = 
arg [ — 14dKb/^d^b ] ' where V is the CKM quark mixing 
matrix [1]. However, the theoretically uncertain contri- 
butions of penguin diagrams (Fig. lb) with different weak 
phases are potentially significant and may shift both the 
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observed CP asymmetries and the branching fractions 
by amounts that depend on the ratios of the penguin to 
tree contributions and their relative phases. A number 
of theoretical estimates exist for the resulting values of 
the branching fractions and CP asymmetries [2~6] . 

The penguin-tree interference in neutral and charged 
B -^ £)(*'£)(*) decays can provide sensitivity to the an- 
gle 7 = arg[-T/,T>;,/\/d\/*i,] [7, 8]. With additional 
information on the branching fractions of -B ^ D)* D^*' 
decays, the weak phase may be extracted, assuming 
SU(3) flavor symmetry between B -^ £)(*)£)(*) and 

B -^ Ds* D^*\ For this analysis, we assume that the 
breaking of SU(3) can be parametrized via the ratios of 
decay constants /^("-j/Zdc*); which are quantities that 
can be determined either with lattice QCD or from ex- 
perimental measurements [9]. 

In addition to presenting measurements of the B'^ -^ 
D(*)+d'-*^- and B+ -^ d(*)+.D(*)o branching fractions, 
and the CP-violating charge asymmetries for the latter 
modes and for B^ -^ D*^D^, we search for the color- 
suppressed decay modes B^ -^ £)(*)^ ]j(*)o ^ which have 
not been previously measured, and determine limits on 
those branching fractions [10]. If observed, the decays B° 
-^ £)(*)" £)(*)o virould provide evidence of VF-exchange or 
annihilation contributions (see Fig. lc,ld). In principle, 
these decays could also provide sensitivity to the CKM 
phase P if sufficient data were available. By combining 
all of these results with information from time-dependent 
CP asymmetries in B^ -^ ]j(*)+]j(*)~ decays and world- 
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FIG. 1: Feynman graphs for B -^ £>(*) [)(*) decays: the tree 
(a) and penguin (b) diagrams are the leading terms for both 
^0 __^ jjM+jjM- and B+ -> _d(*)+^(*)o decays, whereas 
the exchange (c) and annihilation (d) diagrams (the latter of 
which is OZI-suppressed) are the lowest-order terms for B" 
-^ D(*)O0(*''' decays. 



averaged branching fractions of B decays to Dk* D^*^ 
modes, we determine the implications for 7 using the 
method of Refs. [7, 8]. 



II. DETECTOR AND DATA 

The results presented in this paper are based on data 
collected with the BABAR detector [11] at the PEP-II 
asymmetric-energy e+e" collider [12] located at the Stan- 
ford Linear Accelerator Center. The integrated lumi- 
nosity is 210.5 fb^^, corresponding to 231.7 million BB 
pairs, recorded at the T(45') resonance ("on-peak", at a 
center-of-mass (cm.) energy ^/s — 10.58 GeV). 

The asymmetric beam configuration in the laboratory 
frame provides a boost of Pj — 0.56 to the T(4S'). 
Charged particles are detected and their momenta mea- 
sured by the combination of a silicon vertex tracker 



(SVT), consisting of five layers of double-sided detectors, 
and a 40-laycr central drift chamber (DCH), both oper- 
ating in the 1.5-T magnetic field of a solenoid. For tracks 
with transverse momentum greater than 120 MeV/c, the 
DCH provides the primary charged track finding capabil- 
ity. The SVT provides complementary standalone track 
finding for tracks of lower momentum, allowing for re- 
construction of charged tracks with transverse momen- 
tum pt as low as 60 MeV/c, with efficiencies in excess of 
85%. This ability to reconstruct tracks with low pT effi- 
ciently is necessary for reconstruction of the slow charged 
pions from D*+ — > I?°7r+ decays in B ^ £>(*)£>(*) sig- 
nal events. The transverse momentum resolution for the 
combined tracking system is cTpy /pt = 0.0013pT+0.0045, 
where pT is measured in GeV/c. Photons are detected 
and their energies measured by a CsI(Tl) electromag- 
netic calorimeter (EMC). The photon energy resolution 
is as/E = {2.3/£;(GeV)i/'*el.4}%, and their angu- 
lar resolution with respect to the interaction point is 
(Tg ~ (4.2 mrad)/-\/i?( GeV). The measured 7r° mass res- 
olution for TT^'s with laboratory momentum in excess of 

1 GeV/c is approximately 6 MeV/c^. 
Chargcd-particle identification (PID) is provided by an 

internally reflecting ring-imaging Cherenkov light detec- 
tor (DIRC) covering the central region, and the most 
probable energy loss (dE/dx) in the tracking devices. 
The Cherenkov angle resolution of the DIRC is measured 
to be 2.4 mrad, which provides over 5a separation be- 
tween charged kaons and pions at momenta of less than 

2 GeV/c. The dE/dx resolution from the drift chamber is 
typically about 7.5% for pions. Additional information 
to identify and reject electrons and muons is provided 
by the EMC and detectors embedded between the steel 
plates of the magnetic flux return (IFR) . 



III. CANDIDATE RECONSTRUCTION AND B 
MESON SELECTION 

Given the high multiplicity of the final states stud- 
ied, very high combinatorial background levels are ex- 
pected. Selection criteria (described in Sec. Ill A-E) are 
designed to minimize the expected statistical error on 
the B branching fractions (as described in Sec. Ill F). A 
GEANT4-based [13] Monte Carlo (MC) simulation of the 
material composition and the instrumentation response 
of the BABAR detector is used to optimize signal selec- 
tion criteria and evaluate signal detection efficiency. We 
retain sufficient sidebands in the discriminating variables 
to characterize the background in subsequent fits. 



A. Charged track and K" selection 

Charged particle tracks are selected via pattern recog- 
nition algorithms using measurements from the SVT and 
DCH detectors. We additionally require all charged- 
particle tracks (except for those from K*^ -^ tt+tt^ de- 
cays) to originate within 10 cm along the beam axis 



and 1.5 cm in the plane perpendicular to the beam axis 
of the center of the beam crossing region. To ensure 
a wcU-mcasurcd momentum, all charged- particle tracks 
except those from K'^ — > tt+tt" decays and tt"*" from 
£)*+ -^ D^TT^ decays must also be reconstructed from 
at least 12 measurements in the DCH. All tracks that 
meet these criteria are considered as charged pion candi- 
dates. 

Tracks may be identified as kaons based on a likelihood 
selection developed from Cherenkov angle and dE/dx in- 
formation from the DIRC and tracking detectors respec- 
tively. For the typical laboratory momentum spectrum of 
the signal kaons, this selection has an efficiency of about 
85% and a purity of greater than 98%, as determined 
from control samples of D*+ -^ D"7r+, D" — > K^tt^ 
decays. 

We require K^ -^ tt+tt^ candidates to have an invari- 
ant mass within 15 MeV/c^ of the nominal K^ mass [14]. 
The probability that the two daughter tracks originate 
from the same point in space must be greater than 0.1%. 
The transverse flight distance of the K2. from the pri- 
mary event vertex must be both greater than 3cr from 
zero (where a is the measured uncertainty on the trans- 
verse flight length) and also greater than 2 mm. 



B. Photon and tt" selection 

Photons are reconstructed from energy deposits in 
the electromagnetic calorimeter which are not associated 
with a charged track. To reject backgrounds from elec- 
tronics noise, machine background, and hadronic interac- 
tions in the EMC, we require that all photon candidates 
have an energy greater than 30 MeV in the laboratory 
frame and to have a lateral shower shape consistent with 
that of a photon. Neutral pions are reconstructed from 
pairs of photon candidates whose energies in the labo- 
ratory frame sum to more than 200 MeV. The 7r° can- 
didates must have an invariant mass between 115 and 
150 MeV/c^. The 7r° candidates that meet these criteria, 
when combined with other tracks or neutrals to form B 
candidates, are then constrained to originate from their 
expected decay points, and their masses are constrained 
to the nominal value [14]. This procedure improves the 
mass and energy resolution of the parent particles. 



C. Event selection 

We select BB events by applying criteria on the track 
multiplicity and event topology. At least three recon- 
structed tracks, each with transverse momentum greater 
than 100 MeV/c, are required in the laboratory polar 
angle region 0.41 < 6'iab < 2.54. The event must have 
a total measured energy in the laboratory frame greater 
than 4.5 GcV to reject beam-related background. The 
ratio of Fox- Wolfram moments H2/H0 [15] is a param- 
eter between (for "perfectly spherical" events) and 1 
(for "perfectly jet-like" events), and we require this ratio 



to be less than 0.6 for each event, in order to help re- 
ject non-BB background. This criterion rejects between 
30 and 50 percent of non-BB background (depending on 
the decay mode), while keeping almost all of the signal 
decays. 



D. D and D* meson selection 

We reconstruct D^ mesons in the four decay modes 
D° ^ K-TT+, D° -^ K-TT+TT^, DO -^ K-TT+TT-TT+, and 
D° -^ K'^tt^tt^ , and D+ mesons in the two decay modes 
D+ -> K-TT+TT+ and D+ -^ K°tt+ . We require D° 
and D^ candidates to have reconstructed masses within 
±20 MeV/c^ of their nominal masses [14], except for 
£>" -^ K'tt^tt'^, for which we require ±40 MeV/c^ due to 



the poorer resolution for modes containing tt s. These 
criteria correspond to approximately 2.5a of the respec- 
tive mass resolutions. The _D" -^ K~tt^it^ decays must 
also satisfy a criterion on the reconstructed invariant 
masses of the K^t:^ and K^w'^ pairs: the combination of 
reconstructed invariant masses must lie at a point in the 
i^^7r+7r" Dalitz plot [16] for which the expected density 
normalized to the maximum density ( "Dalitz weight" ) is 
at least 6%. Additionally, the daughters of D° and D+ 
candidates must have a probability of originating from a 
common point in space greater than 0.1%, and are then 
constrained both to originate from that common spa- 
tial point and to have their respective nominal invariant 
masses. 

Candidate D*+ and D*° mesons are reconstructed in 
the decay modes D*+ ^ D°n+, D*+ ^ D+tt°, D*" ^ 
D^Tr°, and D*° ^ D"j, using pahs of selected D", D+, 
tt", 7r+, and 7 candidates. The 7r+ from D*+ -^ _D°7r+ 
decays is additionally required to have a cm. momen- 
tum of less than 450 MeV/c. Candidate tt'' mesons from 
D*+ -^ D+tt" and D*° -^ D°7r" are required to have 
cm. momenta p* in the range 70 < p* < 450 MeV/c. 
Photons from I?*° -^ D'^j decays are required to have 
energies in the laboratory frame greater than 100 MeV 
and cm. energies less than 450 MeV. The D* daugh- 
ter particles are constrained to originate from a common 
point in space. After this constraint is applied, the mass 
differences Am of the reconstructed masses of the D* 
and D candidates are required to be within the ranges 
shown in Table I. As shown in Fig. 2, the excellent res- 



TABLE I: Allowed Am{D*- D) ranges for the four D* decay 
modes. 





Minimum 


Maximum 


Mode 


Am (MeV/c^) 


Am (MeV/c") 


D*+ -> D\+ 


139.6 


151.3 


D'+ -. D+^» 


135.0 


146.3 


D*« -^ Z3«^« 


135.0 


149.3 


7?*" ^ D«7 


100.0 


170.0 
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FIG. 2: Distributions of Am in the full data sample for three U* decay modes, 
decays where D'^ decays to K~ ■k'^ . Plot b) shows Am,(D*"'" — D^) for D*"*" 
Plot c) shows Am(L>*" - D°) for D*° -^ D°7r° decays where D'^ decays to K -k 
140.6 MeV/c^, and 142.1 MeV/c^ for the three cases respectively [14]. 



0.142 0.146 
Am (GeV/c2) 

Plot a) shows Am(D*+ -D") for D*+ ^ 

-* D'^Tv^ decays where D^ decays to K'n'^n'^ . 

"^ Nominal values for Am are 145.4 MeV/c^, 



D%^ 



olution in Am for signal candidates makes the Am re- 
quirement a very powerful criterion to reject background 
(see next section), especially for decay modes containing 



E. Variables used for B meson selection 

A S-meson candidate is constructed by combining two 
D^*) candidates that have both passed the selection cri- 
teria described previously. The pairs of D^*^ candidates 
are constrained to originate from the same point in space. 



J 



We form a likchhood variable, i^Mass, that is defined by 
a product of Gaussian distributions for each D mass and 
D* — D mass difference. 

For example, in the decay S" -^ D*'^D*~ , Luass is the 
product of four terms: Gaussian distributions for each D 
mass and double Gaussian (i.e. the sum of two Gaussian 
distributions) terms for each Am term (the D* — D mass 
difference). Defining G(x; fx, a) as a normalized Gaussian 
distribution where x is the independent variable, /i is the 
mean, and a is the resolution, >CMass for B^ -^ D*^ D*~ 
decays is defined as: 



-^Mass = G{mD ; mopoG : '^rno ) X G{m^; m^^^^ , a,n_) X 

/coroG(Amn.+ ; Am^.+ , 0-Ameo,J + (1 - /coro)G(ATO£,.+ ; Am^.+ ,CrAmt,ii) 
/corcG(ATO£,.- ; Amp,.- , CrAm,o„) + (1 - ,/corc)G(ATO£,.- ; ATOr,— , CTAmt^ii) 

PDG PDG 

I 



(2) 



where the subscript "PDG" refers to the nominal 
value [14], and all reconstructed masses and uncertainties 
are determined before mass constraints are applied. For 
(Jmo, we use errors calculated candidate-by-candidate. 
The parameter /core is the ratio of the area of the core 
Gaussian to the total area of the double Gaussian distri- 
bution. This, along with (TAmcoro ^^'^ '^Amta.m is deter- 
mined separately for each of the four D* decay modes 
given above, using MC simulation of signal events that is 
calibrated to inclusive samples of the D* decay modes in 
data. For each of the B decay modes, a higher value of 
'CMass tends to indicate a greater signal likelihood. The 
distributions of — ln(>CMass) for the representative signal 
mode i?° -^ D^D'^ and for the corresponding combinato- 
rial background from generic B'^B^ , B^B^ , cc, and {uu 
+dd +ss) decays, are shown in Fig. 3a. We use Cmuss 



in selecting signal candidates, as will be described in the 
upcoming section. 

We also use the two variables for fully-reconstructed B 
meson selection at the T(45') energy: the beam-energy- 
substituted mass rriES = [(s/2 + Pi ■ Pb)"^ /Ef - p^Y^"^, 
where the initial total e+e^ four-momentum {Ei,pi) 
and the B momentum ps are defined in the laboratory 
frame; and AE = S™ — y/s/2 is the difference between 
the reconstructed B energy in the cm. frame and its 
known value. The normalized distribution of AE for 
the representative signal mode B^ -^ D'^D^ , and for the 
corresponding combinatorial background components, is 
shown in Fig. 3b. 

In addition to >CMass) weSj Eind AE, a Fisher dis- 
criminant T [17] and a £)-meson flight length variable 
L are used to help separate signal from background. The 
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FIG. 3: Distributions of signal selection variables: a) the likelihood variable — ln(/!Mass), b) the /S.E variable, c) the Fisher 
discriminant J-, and d) the D-meson flight length variable L, each for the representative signal mode B" -^ D'^ D^ , and for the 
corresponding combinatorial background from B^B", B^B~ , cc, and (mm +dd +ss) MC simulated decays respectively. In each 
plot, the component distributions are normalized to have the same area below the curves. 



TABLE II: Expected values of the branching fractions B for each B -^ £)(*) £)(*) decay mode, which are used for the purpose of 
determining selection criteria that minimize the expected uncertainty on the measured branching fraction for each mode; also, 
optimized J- and L selection criteria for each mode. An " — " indicates no cut is made in JF or L for that decay mode. 



Mode 



Expected B J^mln Ln 



J— /^v^jc^Lc^i ILJ J mir 

^0 ^ D'+D*- 8.3 X 10"* — — 

B° -^ 0*^0^ 8.8 X 10"* — — 

B" ^ D+D- 3.0 X lO"'' 0.62 1.3 

B° -^ D'^D" 1.0 X 10'^ 0.53 -0.4 



Mode 



■ D*+D'' 



Expected B J^min L„ 



j_jji.jj«jv^ uv^vi /.^ ^ mill -'-'iiiiii 

°B° 1.0 X 10"^ 0.47 -0.4 

i*+_D*'' 1.0 X 10"^ 0.60 — 

B+ -^ D*+_D" 4.4 X 10"" 0.53 -1.3 

B+ -^ B+B*" 4.4 X 10"'' 0.53 0.0 

B+ -^ D+D° 3.0 X 10"'' 0.53 0.5 



B+ 

5 + 



Fisher discriminant assists in the suppression of back- 
ground from continuum events by incorporating infor- 
mation from the topology of the event. The discriminant 
is formed from the momentum flow into nine polar angu- 
lar intervals of 10° centered on the thrust axis of the B 



candidate, the angle of the event thrust axis with respect 
to the beam axis (Ot), and the angle of the B candidate 
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momentum with respect to the beam axis (Ob)' 



j~ — > Cfi^Xi 



(3) 



The values x^ (i = 1, ..., 9) are the scalar sums of the mo- 
menta of all charged tracks and neutral showers in the 
polar angle interval i, xio is |cos0t|, and x\\ is |cos0b|. 
The coefficients a^ are determined from MC simulation 
to maximize the separation between signal and back- 
ground [17]. The normalized distribution of T for the 
representative signal mode _B" -^ D^D^, and for the cor- 
responding background components, is shown in Fig. 3c. 
The flight le ngth var iable L that we consider is defined 
as {ii + l^i)! \/ (y\ + cr|, with the decay lengths li of the 
two D mesons defined as 



^D, = xb + {(.I X pDi) 



(4) 



where xb and xb are the measured decay vertices of the 
D and B, respectively, and po is the momentum of a 
D. The Ui are the measured uncertainties on ii. This 
observable exploits the ability to distinguish the long D 
lifetime. Thus, background events have an L distribution 
centered around zero, while events with real D mesons 
have a distribution favoring positive values. The normal- 
ized distribution of L for the representative signal mode 
S" -^ D'^ D''\ and for the corresponding background com- 
ponents, is shown in Fig. 3d. 



F. Analysis optimization and signal selection 

We combine information from the >CMass, Ai?, !F, 
and L variables to select signal candidates in each 
decay mode. The fractional statistical uncertainty 
on a measured branching fraction is proportional to 
■\/(iV'' + N^)/N^ , where A^" is the number of recon- 
structed signal events and iV'' is the number of back- 
ground events within the selected signal region for a 
mode. The values N^ and N^ are calculated, using de- 
tailed MC simulation of the signal decay modes as well 
as of BB and continuum background decays, by observ- 
ing the number of simulated B decay candidates that 
satisfy the selection criteria for — ln(£Mass), \^E\, T , 
and L. We choose criteria which minimize the expected 
yJiN" + N'')/N'' for each mode. Note that to calculate 
the expected number of signal events N" ^ one must as- 
sume an expected branching fraction, as well as the ratios 
of BB and continuum events using their relative cross- 
sections. These are given, along with the requirements 
on T and i, in Table II. 

For each possible combination of £)*+, D*^ , I?+, and 
D^ decay modes, we determine the combination of selec- 
tion criteria on — ln(>CMass) and | A_B| that minimizes the 
overall expected \/{N'^ + N^)/N^ for each B decay mode 
(see Tables III, IV, and V). The selection criteria for T 
and L are chosen, however, only for each B decay mode 
and not separately for each D^*' mode combination. The 



restrictiveness of the kaon identification selection is also 
optimized separately for each charged and neutral D^*' 
mode. 

Between 1% and 34% of selected B -^ D''*W^*'> events 
have more than one reconstructed B candidate that 
passes all selection criteria in >CMass, ^E, JF, and L, with 
the largest percentages occurring in the decay modes -B" 
-^ £)*°I?*° and 5° -^ D^D'^, and the smallest occurring 
in S" -> D*^D^ and B" ^ D+ D- . In such events, 
we choose the reconstructed B with the largest value of 
'Civiass as the signal B candidate. 



IV. EFFICIENCY AND CROSSFEED 
DETERMINATION 



The efficiencies are determined using fits to rriES dis- 
tributions of signal MC events that pass all selection cri- 
teria in ^Mass, |A-E|, JF, and L. There is a small, but 
non-negligible probability that a signal B decay of mode 
i is reconstructed as a different signal decay mode j. We 
refer to this as crossfeed. Thus, efficiencies can be repre- 
sented as a matrix e^ . where each contributing generated 
event is weighted by the D and D decay mode branching 
fractions. To determine the elements of e^j, we fit the 
mES distributions of signal MC events generated as B 
decay mode i and reconstructed as B decay mode j. The 
distributions are modeled as the sum of signal and back- 
ground probability distribution functions (PDFs), where 
the PDF for the signal is a Gaussian distribution cen- 
tered around the B mass, and the PDF for background 
is an empirical function [18] of the form 



f{x) ex XV 1 — xP' exp[— k(1 — x )], 



(5) 



where we define x = 2mEs/\/s, and k is a parameter 
determined by the fit. In BB MC samples containing 
signal and background decays, we find that the toes dis- 
tribution is well-described by adding a simple Gaussian 
function to the empirical shape in Eq. 5. We fit the toes 
distributions of signal MC events generated as mode i 
and passing selection criteria in mode j to the above dis- 
tribution by minimizing the Xtj of each fit with respect to 
Kij (the K parameter for each mode (i, j)), the number of 
signal events N^; , and the number of background events 
N^j. We determine the efficiencies e^ as NfJN^ , where 
N^ is the total number of signal MC events that were 
generated in mode i. The diagonal elements of the eij 
matrix {i.e. the numbers typically denoted as "efficien- 
cies") are in the range (0.2 - 1.5) x 10^^. The main cross- 
feed source is misidentification between D*^ and D*^ 
candidates. The matrix e^ and the uncertainties on the 
elements of this matrix are given in Table VI. Crossfeed 
between different D submodes {i.e. mode numbers 12-15 
in Table III) is negligible. 
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TABLE III: Key to mode numbers used in Tables IV and V below. 



Mode 


# 


Mode 


# 


D*+ -^ {K-7t+)tt+ 


1 


D'" ^ {K°7r+n-)iT° 


9 


D-+ ^ (K-7T + 7t'>)7T+ 


2 


D*" ^ (X-7r+)7 


10 


D*+ -^ {K-'!t+tt-tt+)tt+ 


3 


D+ -^ K--IT+-IT+ 


11 


D*+ -^ {K°'!t+'!T-)tt+ 


4 


D° ^ K-TV+ 


12 


0"+ -^ (K-7t+tt+)tt° 


5 


D° ^ X-Tr+Tr" 


13 


D'° -^ (X-7r+)7r» 


6 


D« ^ A-7r+7r-7r+ 


14 


D'° ^ (X-7r+7r'')7r'' 


7 


D« ^ X«7r+7r- 


15 


D-° -^ {K-7T+7r-7T + )TV° 


8 







TABLE IV: Optimized — ln(i2Mass) selection criteria used for all B -^ £)(*)£)(*) rnodes. Selected events in a given mode must 
have — ln(/lMass) less than the given value. The D'*' decay modes 1 — 15 are defined in Table III above. Elements with " — " 
above and on the diagonal are modes that are unused since, due to high backgrounds, they do not help to increase signal 
sensitivity. 







D*- 


D*^ 




D 




D' 






12 3 4 5 


6 7 8 9 


10 


11 


12 


13 14 15 




1 


13.0 12.0 17.3 19.8 10.5 


14.6 17.5 9.2 — 


8.9 


8.2 


8.6 


8.5 8.2 8.0 




2 


10.6 11.0 18.3 9.5 


11.5 9.8 10.7 — 


8.7 


8.4 


7.8 


— 8.8 — 


D*+ 


3 


11.7 11.0 9.8 


11.7 9.6 10.4 — 


9.0 


8.8 


9.3 


9.4 9.0 — 




4 








9.6 


15.1 


9.2 












5 


— 


8.2 — — — 


— 


— 


6.6 







6 




12.2 8.4 9.6 7.6 


— 


9.9 


7.6 


6.7 7.2 — 


D*o 


8 
9 
10 




— — 


— 


7.5 
9.2 


5.8 




















D+ 


11 






6.0 


7.3 


5.8 6.5 6.2 




12 










— 


5.2 6.8 — 


D^ 


13 

14 
15 












— 6.2 — 
6.9 — 



TABLE V: Optimized AE selection criteria used for all B — > D^*' O^*' modes. Selected events in a given mode must have \AE\ 
(in MeV) less than the given value. The D*-*' decay modes 1 — 15 are defined in Table III above. Elements with " — " above 
and on the diagonal are modes that are unused since, due to high backgrounds, they do not help to increase signal sensitivity. 







D*- 


D*o 




D- 


DO 






12 3 4 5 


6 7 8 9 


10 


11 


12 13 14 15 




1 


35.5 33.8 30.4 35.2 25.5 


35.7 21.0 26.0 — 


43.6 


18.0 


18.1 20.2 17.1 19.0 




2 


34.5 29.6 23.5 27.4 


40.9 23.9 21.4 — 


29.3 


19.4 


25.9 — 19.5 — 


JJ*+ 


3 

4 


23.5 23.7 18.2 


34.0 30.6 20.6 — 


27.3 


18.6 
21.9 


19.0 20.4 17.1 — 
16.9 19.7 — — 












5 


— 


19.1 — — — 


— 


— 


16.4 — — — 




6 

7 
8 




35.1 23.0 27.3 25.5 


— 


23.9 
20.0 
16.6 


17.4 19.6 17.4 — 


D*^ 















9 




— 


— 


— 


24.5 — — — 




10 






— 


— 


— — — — 


D+ 


11 






15.1 


15.5 19.2 15.4 15.5 




12 










— 18.7 16.1 — 


DO 


13 










— 19.0 — 


14 










15.9 — 




15 










— 
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TABLE VI: Elements of the efficiency and crossfeed matrix e^j, and their respective uncertainties, used to calculate the 
branching fractions and charge asymmetries, as described in the text. All values are in units of 10""*. Uncertainties on the 
last digit(s) are given in parentheses. Elements with " — " correspond to values that are zero (to three digits after the decimal 
point). The column corresponds to the generated mode and the row corresponds to the reconstructed mode. 



Mode 



D'+D*- D'-^D'f D+D- D*°D*° D*°D° 73°D" D*+ D*° D*+ D"^ D+ D*" D+D 



+ no 



D*+D*- 
D+D- 

^.000 

D*+W 
D+D*° 
D+W 



14.24(6) 0.010(3) — — — — 0.18(1) _ _ _ 

0.020(3) 11.52(6) — — — — 0.010(3) 0.040(3) 0.08(1) — 

— — 9.51(8) ______ 0.010(3) 

0.080(3) — — 2.60(2) 0.030(3) — 0.42(1) 0.010(3) — — 

— — — 0.020(3) 3.40(2) — 0.010(3) 0.46(1) 0.010(3) — 

— — — — 0.010(3) 12.02(10) — 0.010(3) 0.020(3) — 
2.60(2) — — 0.23(1) 0.010(3) — 7.52(4) 0.07(1) — — 

0.040(3) 0.06(2) — — 0.11(5) — 0.03(2) 13.51(25) 0.040(3) — 



0.41(1) 
0.020(3) 0.06(1) 



0.010(3) 0.010(3) 



0.050(3) 



0.070(3) 3.70(3) — 
0.010(3) 0.020(3) 14.93(9) 



V. BRANCHING FRACTION RESULTS 

In order to determine the number of signal events in 
each mode, one must not only account for background 
which is distributed according to combinatorial phase 
space, but also for background which can have a differ- 
ent distribution in mES- It is possible for a component 
of the background to have an tties distribution with a 
PDF that is more similar to signal {i.e. a Gaussian dis- 
tribution centered around the B mass) than to a phase- 
space distribution. Such a component is known as "peak- 
ing" background and typically derives from background 
events that have the same or similar final state particles 
as the signal decay mode. For example, in B^ -^ D^D~ , 
peaking background primarily comes from the decays 
B" ^ DKX or 5° -^ DnX, where D -^ Rtttt and 
X is K'^, p, ai or cu, and the light mesons (KX) or {ttX) 
fake a, D —^ Kmr decay. The optimization procedure 
that was detailed in Sec. Ill F eliminates decay submodes 
that have a large enough amount of peaking (in addition 
to combinatorial) background to decrease, rather than 
increase, the sensitivity for a particular decay; the final 
selection was detailed in Tables II, IV, and V. We de- 
termine the amount of peaking background Pi in each B 
decay mode i via fitting the toes distributions of BB MC 
simulated events. We minimize the xf of each fit, allow- 
ing the variables nf (representing the "ARGUS param- 
eter" described earlier) , the number of expected peaking 
background events in data Pi, and the number of phase- 
space background events N- ^, to float. The fitted 
number of peaking background events Pi is compatible 
with zero, within two standard deviations, for all modes 



similar to the PDFs used for efhcicncy and peaking back- 
ground fits described above. We minimize the xf of each 
data fit, allowing the parameter Ki, the number of signal 
events in data N^^^, and the number of background events 
in data N^ *, each to float. The toes distributions and 
the results of the fits are shown in Fig. 4. The branching 
fractions Bi are then determined via the equation 



j 



EijBjNs 



n! 



P 



(6) 



where Nb = N^g = (231. 7± 2.6) x 10'^ is the total num- 
ber of charged or neutral B decays in the data sample, 
assuming equal production rates of charged and neutral 
B pairs. 

We determine the branching fractions as 



B.^J2'^(Nf'-P^)/N' 



(7) 



(where e~ is the inverse of matrix e^ ) yields the branch- 
ing fractions given in Table VII. Note that the measured 
branching fractions for the three modes B^ -^ £)(*)0£)(*)o 
are not significantly greater than zero. Thus, we have 
determined upper limits on the branching fractions for 
these modes. The 90% confidence level (C.L.) upper 
limits quoted in Table VII are determined using the 
Fcldman- Cousins method [19] and include all systematic 
uncertainties detailed below. Since the branching frac- 
tions can be correlated through the use of Eq. 6, we also 
provide the covariance matrix, with all systematic uncer- 
tainties included, in Table VIII. The covariance matrix 
is obtained via the approximation given in [20] . 



We then fit the actual data to determine the number 
of reconstructed signal events in each mode. We fit the 
to-es distributions of reconstructed B decays that pass 
all selection criteria in each mode i to a sum of a Gaus- 
sian distribution and a phase space distribution (Eq. 5), 



VI. BRANCHING FRACTION SYSTEMATIC 
UNCERTAINTIES 

Table IX shows the results of our evaluation of the 
systematic uncertainties on the branching fraction mea- 
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FIG. 4: Distributions of tues for selected candidates in each D**'_D'*-' mode. The error bars represent the statistical errors 
only. The solid lines represent the fits to the data, and the shaded areas the fitted background. 
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TABLE VII: Results of the fits for the ten signal decay modes: the number of events for fitted signal A'''"^, the peaking back- 
ground P, and the crossfeed C, the branching fractions B, 90% C.L. upper limits on branching fractions, previous measurements 
of branching fractions (for modes that have previous measurements), and charge asymmetries. The uncertainties are statistical. 
For the final branching fraction and charge asymmetry results, the systematic errors are also given. 



Mode 




JV^is 


P 


C 


B (10-^) 


U.L. 

(10-^) 


Previous B 
results (10"") 


Acp 


B° -* 


D'+D*- 


270±19 


-1±2 


4±1 


8.1±0.6±1.0 




8.1 ±0.8 ±1.1 [21] 

8.3 ±1.6 ±1.2 [22] 

9.9+^.^3 ± 1-2 [23] 




B° ^ 


D'^D'f 


156±17 


1±3 


2±1 


5.7±0.7±0.7 




8.8 ±1.0 ±1.3 [24] 
11.7±2.6l|^ [25] 
6.7ti;? ± 1.1 [26] 


0.03 ± 0.10 ± 0.02 [27] 


B° ^ 


D+D- 


63±9 


1±2 


0±0 


2.8±0.4±0.5 




1.91 ±0.51 ±0.30 [28] 




B" ^ 


D*°D'° 


0±6 


-2±2 


0±0 


-1.3±1.1±0.4 


0.9 


< 270 [29] 




B" -. 


D*°D" 


10±8 


-2±3 


1±1 


1.0±1.1±0.4 


2.9 






B" ^ 


D°D° 


-11±12 


-8±4 


0±0 


-0.1±0.5±0.2 


0.6 






B+-> 


jj* + ^*0 


185±20 


-5±4 34±4 


8.1±1.2±1.2 




W.5tll±2.0 [26] 
< 110 [29] 


-0.15 ± 0.11 ±0.02 


B+ ^ 


D*+D° 


115±16 


1±4 


3±1 


3.6±0.5±0.4 




4.57 ±0.71 ±0.56 [28] 
< 130 [29, 30] 


-0.06 ± 0.13 ±0.02 


B+ ^ 


D+D*° 


63±11 


3±3 


9±2 


6.3±1.4±1.0 




< 130 [29, 30] 


0.13 ± 0.18 ± 0.04 


B+ ^ 


D+D° 


129±20 


-2±5 


1±1 


3.8±0.6±0.5 




4.83 ± 0.78 ± 0.58 [28] 
< 67 [29] 


-0.13 ±0.14 ±0.02 



TABLE VIII: Covariances of B — > D^*'D^*' branching fractions (with all systematic uncertainties included), in units of 10 



Mode D*+D*- 


D*±D^ 


D+D- 


^.00*0 


^.0^0 


D°D" 


D*+_D*« 


D*+^o 


D+D*'> 


D+D° 


D*+D*- 1.26 


0.55 


0.22 


-0.15 


0.07 


-0.01 


0.73 


0.33 


0.54 


0.30 


D*±D^ 


0.91 


0.26 


-0.08 


0.04 


-0.01 


0.46 


0.19 


0.37 


0.26 


D+D- 




0.39 


-0.03 


0.02 


0.00 


0.16 


0.08 


0.26 


0.16 


D*°D'-° 






1.27 


-0.04 


0.00 


-0.53 


-0.06 


-0.13 


-0.05 


D*°D° 








1.25 


0.00 


0.07 


-0.02 


0.05 


0.02 


D°D° 










0.22 


-0.01 


0.00 


-0.01 


0.00 


D*+D"' 












2.60 


0.31 


0.55 


0.27 


D*+W 














0.43 


0.19 


0.11 


D+D*° 
















2.61 


0.27 


D+D" 


















0.53 



surcmcnts. 

a. Submode branching fractions The central values 
and uncertainties on the branching fractions of the D 
and D* mesons are propagated into the calculation of 
the branching fraction measurements. The world average 
measurements [14] are used. 

b. Charged track finding efficiency From studies of 
absolute tracking efficiency, we assign a systematic uncer- 
tainty of 0.8% per charged track on the efficiency of find- 
ing tracks other than slow pions from charged D* decays 
and daughters of K'^ decays. For the slow pions, we as- 



sign a systematic uncertainty of 2.2% each, as determined 
from a separate efficiency study (using extrapolation of 
slow tracks found in the SVT into the DCH tracking 
detector and vice- versa). Track finding efficiency uncer- 
tainties are treated as 100% correlated among the tracks 
in a candidate. These uncertainties are weighted by the 
D and D* branching fractions. 

c. K'^ reconstruction efficiency From a study of the 
K'^ reconstruction efficiency (using an inclusive data 
sample of events containing one or more K'^, as well as 
corresponding MC samples), we assign a 2.5% systematic 



16 



TABLE IX: Estimates of branching fraction systematic uncertainties (as percentages of the absolute values of the branching 
fraction central values) for all B modes, after propagating the errors through Eq. 6. The totals are the sums in quadrature 
of the uncertainties in each column. Note that the term "Dalitz weight" refers to the selection on the reconstructed invariant 
masses of the K^tv^ and if^Tr" pairs for D" -^ A'^tt+tt" decays that was described in Sec. HID. 



Mode 


D*+D*- 


D*±Z)T 


D+D- 


D*o_D*o 


D*''W 


D^'W 


^.+ 0*0 


D'+W 


D+^*o 


D+W 


D*+ BFs 


1.4 


0.7 


0.0 


0.9 


0.1 


0.0 


0.7 


0.7 


0.0 


0.0 


D*° BFs 


0.0 


0.0 


0.0 


4.9 


1.6 


0.0 


2.1 


0.0 


4.4 


0.0 


D" BFs 


5.0 


2.7 


0.0 


7.4 


3.7 


5.7 


5.2 


4.5 


3.3 


2.7 


D+ BFs 


1.4 


6.5 


13.2 


0.1 


0.2 


0.4 


0.1 


0.3 


6.5 


6.5 


Tracking efficiency 


7.9 


6.5 


4.8 


7.9 


3.0 


4.7 


6.0 


6.0 


3.8 


4.4 


Kg efficiency 


0.3 


0.2 


0.0 


0.0 


0.1 


0.0 


0.1 


0.2 


0.3 


0.2 


Neutrals efficiency 


2.5 


1.0 


0.0 


8.4 


2.9 


1.9 


4.6 


1.6 


4.3 


1.0 


Kaon identification 


4.6 


4.7 


5.0 


7.3 


4.9 


5.4 


5.0 


4.6 


4.6 


4.7 


-CMass cut 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


JT cut 


0.0 


0.0 


0.9 


0.9 


0.9 


0.9 


0.9 


0.9 


0.9 


0.9 


L cut 


0.0 


0.0 


0.8 


0.8 


0.8 


0.8 


0.0 


0.8 


0.8 


0.8 


AE cut 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


Dalitz weight cut 


1.0 


0.5 


0.0 


1.4 


0.2 


1.0 


1.0 


0.8 


0.7 


0.5 


P(X') cut 


3.8 


3.8 


3.8 


3.8 


3.8 


3.8 


3.8 


3.8 


3.8 


3.8 


Fit model 


1.8 


3.6 


3.1 


5.4 


6.7 


44.6 


4.9 


2.8 


7.0 


3.6 


Spin alignment 


1.0 


0.0 


0.0 


6.1 


0.0 


0.1 


4.1 


0.0 


0.0 


0.0 


Peaking background 


0.9 


2.0 


2.9 


24.5 


32.3 


144.6 


3.1 


3.4 


4.9 


4.0 


Crossfeed 


0.4 


0.6 


0.8 


1.9 


1.1 


1.6 


0.6 


0.4 


1.0 


0.6 


N.J, 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


Total 


12.0 


12.3 


16.1 


31.0 


34.2 


151.7 


13.6 


11.0 


14.8 


11.9 



uncertainty for all modes containing a K^,- The value 
2.5% comes from the statistical uncertainty in the ratio 
of data to MC yields and the variation of this ratio over 
different selection criteria. The uncertainty is weighted 
by the D and D* branching fractions. 

d. tt'^ and 7 finding efficiency From studies of the 
neutral particle finding efficiency through the ratios of 
T+ -^ p^ {tt'^ TT^)^ to r"*" ^- TT+i^ between data and MC, 
we assign a 3% systematic uncertainty per tt^, includ- 
ing the slow 7r° from D* and £>*" decays. For isolated 
photons from D*" decays, we assign a 1.8% systematic 
uncertainty, 100% correlated with the 7r° efficiency un- 
certainty. These uncertainties are weighted by the D and 
D* branching fractions. 

e. Charged kaon identification We assign a system- 
atic uncertainty of 2.5% per charged kaon, according to 
a study of kaon particle identification efficiency (using 
kinematically-reconstructed D^ -^ K^n^ candidates). 
The uncertainty is weighted by the D and D* branching 
fractions. 

/. Other selection differences between data and MC 
Differences in momentum measurement, decay vertex 
finding efficiency, etc., can result in additional differences 
between efficiencies in data and in MC. We use a sample 
of the more abundant i?° -^ D*~^D*~ events in data, se- 
lected in a similar manner as the B -^ D^*'D^*' modes, 
to determine these uncertainties. To estimate the sys- 



tematic error arising from differences between the data 
and MC D and D* mass resolutions, we calculate the 
number of D^D* events seen in the data and MC as a 
function of the Z^Mass cut, while fixing the other selection 
criteria to their nominal values. The number of observed 
events is extracted from a fit to the tties distribution. We 
then plot the ratio of the data yield (A'^data) to the MC 
yield (Nmc) as a function of the £Mass cut over a range of 
values that gives the same efficiencies as in the £)(*)£)(*) 
analyses. We find the rms of the Ndata/^MC ratio and 
assign this as a systematic uncertainty for applying this 
cut. The same technique is used to determine the sys- 
tematic uncertainties from all other selection criteria in 
Table IX: the selections on !F, L, AE, the reconstructed 



("Dalitz weight"). 



invariant masses for Z?" -^ K tt - 
and vertex P(x^). 

g. Fit model The data yield is obtained from an 
TTiES fit where the mean (n) and width {a) of the B mass 
and the end-point (y^/2) of the phase-space distribution 
(Eq. 5) are fixed. These parameters are estimated and 
have associated uncertainties. The nominal value of a is 
determined from signal MC for each B decay mode. To 
estimate the systematic uncertainty due to possible dif- 
ferences between the tties resolutions in data and signal 
MC, we first look at this difference (Act ~ crdata — ctmc) 
for those modes with high purity, including our control 
sample. These differences are consistent with zero, jus- 
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tifying our use of (Tmc in obtaining the data yield. We 
then find the weighted average of Atr, which is given by 
(0.11 ± 0.08) MeV/c^. As a conservative estimate, we re- 
peat the data yield determinations by moving a up and 
down by 0.2 MeV/c^, and take the average of the ab- 
solute values of the changes in each data yield as the 
systematic uncertainty of fixing a to the MC value for 
that B mode. A combined fit of common modes in data 
is used to determine the nominal values for \i and for the 
endpoint of the rriES distribution \/s/2. Hence, we move 
the parameters up and down by their fitted errors (0.2 
MeV/c^ for ^ and 0.1 MeV/c^ for v^/2) to obtain their 
corresponding systematic uncertainties. The quadratic 
sum of the three uncertainties from /i, a and \/s/2 gives 
the systematic uncertainty of the fit model for each B 
mode. 

h. Spin- alignment dependence The B^ -^ D*^D*^ , 
B° -^ D*°D*°, and B+ -^ d*°D*+ decays are pseu- 
doscalar — > vector vector (VV) transitions described 
by three independent helicity amplitudes Aq, Ah, and 
A± [31]. The lack of knowledge of the true helicity am- 
plitudes in the B -^ VV final states contributes a system- 
atic uncertainty to the efficiency. The dominant source of 
this effect originates from the p-r-dependcnt inefficiency 
in reconstructing the low-momentum "soft" pions in the 
D*"*" and D*" decays, and the fact that the three helicity 
amplitudes contribute very differently to the slow pion 
Pt distributions. To estimate the size of this effect, MC 
samples are produced with a phase-space angular distri- 
bution model for the decay products. Each event is then 
weighted by the angular distribution for given input val- 
ues of the helicity amplitudes and phase differences. The 
efficiency is then determined for a large number of ampli- 
tude sets and the observed distributions in efficiencies are 
used to estimate a systematic uncertainty. For a given 
iteration, a random number, based on a uniform PDF, is 
generated for each of the three parameters: R±,a, and 
77, where 



Ra 



l^op-|A||P 



l^ol 



l^llp 



\Aj 



I^IlP' 



(8) 



and 77 is the strong phase difference between ^0 and A\\. 
Since R± for B'^ -^ D*^D*^ has already been mea- 
sured [36], a Gaussian PDF with mean and width fixed 
to the measured values is used instead for that mode. 
The events of the MC sample are weighted by the cor- 
responding angular distribution and the efficiency is de- 
termined (after applying all selection cuts) by fitting the 
rriES distribution and dividing by the number of gener- 
ated events. The procedure is repeated 1000 times for 
each B -^ VV sample. The relative spread in efficien- 
cies (rms divided by the mean) is used to estimate the 
systematic uncertainty due to a lack of knowledge of the 
true amplitudes. 

i. Peaking background and crossfeed The uncertain- 
ties on the peaking background vector Pi and on the 
efficiency matrix e^ are dominated by the available MC 
statistics. The resulting uncertainties on each clement 



of the vector and matrix are propagated through to the 
branching fraction results via the formalism of Eq. 6. 

j. Number of BB The number of BB events in the 
full data sample, and the uncertainty on this number, 
arc determined via a dedicated analysis of charged track 
multiplicity and event shape [15]. The uncertainty in- 
troduces a systematic uncertainty of 1.1% on each of the 
branching fractions. 



VII. MEASUREMENT OF CP- VIOLATING 
CHARGE ASYMMETRIES 

To obtain the charge asymmetries Acp (defined in 
Eq. 1), we perform unbinned extended maximum like- 
lihood fits to the Toes distributions of the selected events 
in each of the four charged-_B decay modes D*+I?*°, 
£)*+£)°, I?+D*°, D+D°, and their respective charge con- 
jugates, and in the neutral- -B decay mode D* D^ , using 
Eq. 5 as the PDF for the combinatorial background for 
both charges in each pair. The free parameters of each of 
the five fits individually are: 1) the combinatorial back- 
ground shape parameter k, 2) the total number of signal 
events, 3) the total number of background events, and 
4) the "raw" charge asymmetry A. Parameters 1 and 
3 are considered (and thus constrained to be) the same 
for both charge states in each mode; this assumption is 
validated in MC simulation of the background as well as 
in control samples of B° -^ D*~ p'^ and B^ -^ D*~af 
decays in data. The results of the fits are shown in Fig. 5. 
Two potentially biasing effects must be considered: there 
can be a asymmetry in the efficiencies for reconstructing 
positively- and negatively-charged tracks, and peaking 
background and crossfeed between the modes can cause 
a small difference between the measured ("raw") asym- 
metry and the true asymmetry. The former of those two 
effects is discussed in Sec. VIII below. Regarding the 
latter, to obtain the charge asymmetries Acp from the 
"raw" asymmetries A, very small corrections for peak- 
ing background and crossfeed between modes must be 
made. Using the terminology of Eq. 6, and considering 
the branching fractions Bi to be sums of a "+" mode 
(with a i?° or i?"*", containing a b quark, as the initial 
state) and a "— " mode (with a B" or B~ , which contain 
a b quark, as the initial state): Bi = Bf + B^ , we have 
the two equations 



J2e,,BfNB^Nf'^-Pt 



(9) 



for the "+" and "— " states respectively, which imply 

^ e^mr- Pi) ± (<^+ - P;)]/Nb. (10) 



As 



r,r - r+ B7 - b7 



^CPA 



rr + r+ Br + Bf 



(11) 
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5.2 5.22 5.24 5.26 5.28 5.3 5.2 5.22 5.24 5.26 5.28 5.3 



m^3(GeV/c') 



m^g(GeV/c') 



FIG. 5: Fitted distributions of rriES for the two conjugate states of each of the five relevant modes. The error bars represent 
the statistical errors only. The solid lines represent the fits to the data, and the shaded areas the fitted background. The raw 
asymmetries A are the normalized differences in the amount of signal between the members of each conjugate pair. 
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we have 



Eie-'mf'-pn-iK'^-pn 






Since Nf^ = Nf^' 



T..^'mr--P3 



(N. 



sig+ 



pn 



in a mode Aj 



w. 



Ar"K-+Ar' 



A^^'^ and the "raw" asymmetry 
ijqr, we have 



AcP.i ~ 



E 



3 ^V 



'\A,Nf<^- 



AfPj] 



E 



3 ^ij 



[iVfs 



^.] 



(13) 



where Aj 



p^-p+ 

-^ — -^ are the charge asymmetries of the 

j j 

peaking backgrounds. The total yields iV^'®, peaking 
backgrounds Pj, and efhciency matrix Cij are identical 
to those used for the branching fraction measurements 
and are given in Tables VII and VI. The values A^ are 
nominally set to and are varied to obtain systematic 
uncertainties due to the uncertainty on the charge asym- 
metry of the peaking background (see Sec. VIII). Thus, 
Eq. 13 is used to determine the final Acp values from 
the measured asymmetries, in order to account for the 
small effects due to peaking background and crossfeed 
between modes. The measured Acp values are given in 
Table VII. They are all consistent with zero, and their 
errors are dominated by statistical uncertainty. 



uncertainty of 0.5% for each Acp measurement for modes 
with a charged slow pion. 

h. Charge asymmetry from tracks other than slow it 
Auxilliary track reconstruction studies place a stringent 
bound on detector charge asymmetry effects at transverse 
momenta above 200 MeV/c. Such tracking and PID sys- 
tematic effects were studied in detail in the analysis of 
B -^ (t)K* [34]. We assign a 0.2% systematic per charged 
track, thus an overall systematic of 0.4% per mode (as 
the positively charged and negatively charged decays for 
each mode have, on balance, one positive vs. one nega- 
tive track respectively). This systematic uncertainty is 
added linearly to the slow pion charge asymmetry sys- 
tematic due to potential correlation. 

c. Amount of peaking background Peaking back- 
ground can potentially bias Acp measurements in two 
ways: 1) a difference in the total amount of peaking back- 
ground from the expected total amount can, to second 
order, alter the measured asymmetry between the posi- 
tively charged and negatively charged decays, 2) a more 
direct way for peaking background to alter the measured 
Acp would be if the peaking background itself were to 
have an asymmetry between the amount that is recon- 
structed as positively charged and the amount recon- 
structed as negative. 1) is discussed here; 2) is discussed 
in the paragraph below. The systematic uncertainty due 
to the uncertainty on the total amount of peaking back- 
ground in the five decays is determined via the formalism 
of Eq. 13. Namely, the uncertainty is given by 



VIII. SYSTEMATIC UNCERTAINTIES ON 
CHARGE ASYMMETRY MEASUREMENTS 

Table X shows the results of our evaluation of the vari- 
ous sources of systematic uncertainty that are important 
for the Acp measurements. 

a. Slow TT charge asymmetry A charge asymme- 
try in the reconstruction efhciency of the low-transverse- 
momentum charged pions from D -^ D^tt^ decays can 
cause a shift in Acp by biasing the rates of positively 
charged vs. negatively charged decays for each mode. 
We estimate this systematic uncertainty by using data 
control samples of S° -> D*' X+ and B° -^ D*+X- 
decays, where X is either tt, p, or ai, and determining if 
there is an asymmetry in the number of Z?*+ vs. D*^ 
reconstructed. There are two potential biases of this 
technique: 1) a charge asymmetry in tracks other than 
the slow charged pions, and 2) the presence of doubly- 
Cabibbo-suppressed B^{B'^) -^ D*^X^ decays which 
could potentially introduce a direct-CP-violating asym- 
metry between the two states in the control sample. Dis- 
cussion of 1) is detailed in the paragraph below, and 
the rate of 2) has been determined in analyses such as 
Refs. [32] and [33] to be of order 0.1%, well below the 
sensitivity for this measurement. We combine the infor- 
mation from the control sample modes and determine an 



5A, 



{E,er^'A,Nf^)xJj:^{e-'y{SP) 



^3^"] 



CP,i 



(E,^.^-M<^-^,])^ 



(14) 



where {6P)j are the uncertainties on the amount of peak- 
ing background (which are given, along with the other 
parameters in the equation, in Table VII). 

d. Acp of peaking background The systematic un- 
certainty due to the Acp of the peaking background is 
also determined using the formalism of Eq. 13. Namely, 
the uncertainty is given by 



S'A. 



CPA 



CL,ellA,Nf-)^JY.,{^^?{5APy]Pj 



(E 



3 V ^3 



P3])iE,e^'A,N[ 



(15) 



Investigation of the sources of the peaking background in 
these modes motivates a conservative choice of 0.68 for 
the {6A^)j values. 

e. Amount of crossfeed The systematic error due 
to uncertainties in the amount of crossfeed between the 
modes is also determined via the formalism of Eq. 13. 
Namely, the uncertainty is given by 



TABLE X: Summary of the systematic uncertainties estimated for the Acp asymmetries, in %. 
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Systematics source 



B» 



D*±DT B^ 



D*+D*o B^ 



D'+D° B+ -^ D+D*° B^ 



D+D° 



Slow pion charge asymmetry 


0.53 


0.53 


0.53 


— 


— 


Charge asymmetry from other tracks 


0.80 


0.80 


0.80 


0.80 


0.80 


Amount of peaking bkgd. 


0.06 


0.42 


0.19 


0.64 


0.53 


Acp of peaking bkgd. 


0.42 


0.09 


0.58 


3.36 


0.85 


Crossfeed uncertainty 


0.00 


0.00 


0.01 


0.31 


0.00 


rriES resolution uncertainty 


0.20 


0.04 


0.04 


0.01 


0.14 


B mass uncertainty 


0.20 


0.37 


1.38 


1.38 


0.53 


Uncertainty in y'i 


0.00 


0.03 


0.08 


0.05 


0.05 


Potential fit bias 


0.74 


f.97 


1.19 


0.53 


1.66 



TOTAL 5{Acp) 



1.6 



2.4 



2.3 



3.8 



2.2 



S"A. 



CP,i 



E,-,A-A^fcov(e--\e-,^)AiVr 



(i:,^^^-4,<^)x(v/E,.[< 



-P,]cov(er.i,,-i)[ivf-Pfe] 



(E 



e--M<" 



] y 



PA? 



(16) 



r 



The covariance between the elements of the inverse effi- 
ciency matrix is obtained using the method of Ref. [20] . 
The very small systematic uncertainty due to crossfeed is 
thus obtained using Eq. 16 and the amounts of crossfeed 
and their uncertainties that are given in Table VI. 

/. Uncertainty in tties resolution, B mass, and ^/s 
The uncertainties in toes resolution and the beam en- 
ergy a/s are determined by varying these parameters 
within their fitted ±lcr ranges and observing the resulting 
changes in Acp- The uncertainty in the reconstructed B 
mass can also have an impact on the fitted ttt-es distribu- 
tions and thus on the fitted Acp values. Varying the B 
mass between the fitted value and the ilcr range of the 
nominal i?° or B^ invariant mass allows the determina- 
tion of the resulting effect on the Acp values. 

g. Potential fit bias Uncertainties in the potential 
biases of the Acp fits are determined by performing the 
fits on large samples of MC simulation of the signal de- 
cay modes and of BB and continuum background decays. 
All results are consistent with zero bias, and the uncer- 
tainties of the fitted asymmetries on the simulated data 
samples are conservatively assigned as systematic uncer- 
tainties from biases of the fits. 



IX. IMPLICATIONS FOR 7 

Information on the weak phase 7 may be obtained by 
combining information from B -^ £){*)£)(*) and B -^ 



jjyi £){*) branching fractions, along with CP asymmetry 
measurements in i? — > D'-*^D^*\ and using an SU(3) re- 
lation between the D^*'^D^*^ and d'i'^D^*^ decays [7, 8]. 
For this analysis, we assume that the breaking of SU(3) 
can be parametrized via the ratios of decay constants 
/r)(*)//u(*) J which are quantities that can be determined 
either with lattice QCD or from experimental measure- 
ments [9]. 

In this model, one obtains the relation (for i?'^ -^ 
D+D- and individual hclicity states of S" ^ D*+D*-): 



A^ 



where 



Qfl cos(2/3 + 27) - ttindir sin(2/3 + 27) - g 
cos 27 — 1 



B^-^i\A 



D\2 



ttdii 



^indii 



A^ + A^ 



f 2ActAut cos S cos 7, 



-2^ct-4„tsin(5sin7, 



^e~''P{A'^)*A^] = -Aitsm2p 
—2ActAut cos (5 sin (2/3 + 7) 
-^2^sin(2/3 + 27), 



and 



2 



B' 



(17) 



(18) 



(19) 



(20) 



(21) 
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A^ and A^ represent amplitudes of a given B'^ and 
B^ -^ £)(*)+£)(*)- decay respectively, B represents the 
corresponding average branching fraction, and Odir and 
flindir represent the corresponding direct and indirect 
CP asymmetries respectively. The phases /3 and 7 are 
the CKM phases and 6 is a strong phase difference. 
Act = \iT + E + Pc- Pt- PEwWclVcdl and Ant = 
\{Pu -Pt- PEwW-i^b^'^dl are the magnitudes of the com- 
bined B ^> D^*^D'-*^ decay amplitudes containing V*ijVcd 
and V*^ud terms respectively, and the T, F, and E 
terms are the tree, penguin, and the sum of exchange 
and annihilation amplitudes respectively [7] . One can di- 
rectly measure the parameters ,B, adir, and aindir using in- 
formation from B -^ £)(*)£)(*) decays; the parameter Act 
using information from B -^ Ds* D^*^ decays; and the 
weak phase /? can be obtained from the measurements of 
sin 2/3 based on _B" -^ ccK^ decays [35] thus allowing for 
solution of 7 (up to two discrete ambiguities) via Eq. 17. 
As the vector-pseudoscalar modes _B" -^ D*^D^ are not 
CP eigenstates, a slightly more complicated analogue to 
Eq. 17 is needed for these modes [8]. Measurement of 
Acp for D*^D^ is also necessary to obtain information 
on 7 from the vector-pseudoscalar modes. 

Using these relations, there are four variables besides 
(3 for each B -^ £)(*)!)(*) decay for which to solve: Act, 
Auti <5, and 7. The branching fraction and the direct and 
indirect CP asymmetries of the B -^ £)(*)£)(*) decay pro- 
vide three measured quantities. The other measurement 
that can be used is the branching fraction of the cor- 
responding B — > Ds D^*> decay, by using the relation 
expressed in Eq. 22. 

The values aindir can, of course, only be measured in 
the neutral B — > £)(*)£)(*) decays. However, the charged 
B -^ £)(*)£)(*) decays can supplement the neutral de- 
cays by adding information on B and adir, assuming 
only isospin symmetry between the charged and neutral 
modes. Thus, information from the charged B decay 
modes can assist the 7 determination. 

SU(3)-breaking effects can distort the relation between 

i:)(*)Z)(*) and Di*''D'-*^ decays as expressed in Eq. 17. 
However, the SU(3)-breaking can be parametrized by the 
ratio of decay constants f jji*) / f d^') , such that the am- 
plitude for B -f Di*''D'-*^ decays 



-^'ct = fni*)lfD(-) X At/ sin 0c 



(22) 



where Oc is the Cabibbo angle [14] and the parentheses 
around the asterisks correspond to the B -^ £)(*)£)(*) 
and B -^ Dg D'^*'> decays that are used. The theoretical 
uncertainty of this relation is determined to be 10% [7]. 

We thus use the information from the vector-vector 
(VV) decays S" ^ D*+D*- and B+ -> D*+D*° and 
pseudoscalar-pseudoscalar (PP) decays B^ -^ D^D^ 
and 5+ -^ D^D'^, as well as the vector-pseudoscalar 
(VPj_ decays B° ^ D*^D^, B+ -^ D*+W, and B+ -^ 
D^ D*^ , to form constraints on 7 using the method of 
Refs. [7, 8]. 

To use the VV decays, we must make the assump- 



tion that the strong phases for the and || helicity am- 
plitudes are equal. The constraints from the PP de- 
cays require no such assumption. The assumption of 
equal and || helicity amplitudes is theoretically sup- 
ported by a QCD factorization argument described in [8]. 
Then, using Eq. 17, we combine the B'^ -^ D*^D*^ 
and 5+ — > D*^D*^ branching fractions and Acp in- 
formation given above with measurements of the B^ 
-^ D*-D*+ and B+ -^ D*+D*° branching fractions [14], 
measurements of the B^ -^ D*^D*~ time-dependent CP 
asymmetries [21, 36], and the world-average values of 
sin 2/3 [35] and sin 6*^ [14]. 

We use a fast parametrized MC method, described in 
Ref. [8], to determine the confidence intervals for 7. We 
consider 500 values for 7, evenly spaced between and 
27r. For each value of 7 considered, we generate 25000 
MC experiments, with inputs that are generated accord- 
ing to Gaussian distributions with widths equal to the 
experimental errors of each quantity. For each exper- 
iment, we generate random values of each of the ex- 
perimental inputs according to Gaussian distributions, 
with means and sigmas according to the measured cen- 
tral value and total errors on each experimental quan- 
tity. We make the assumption that the ratio /d»//d* 
is equal to fajfo = 1-20 ± 0.06 ± 0.06 [9], allowing for 
the additional 10% theoretical uncertainty [7]. We then 
calculate the resulting values of Act, Odir, Oindir, and B, 
given the generated random values (based on the exper- 
imental values). When the quantities Cdir, Oindir, and B, 
along with (3 and the value of 7 that is being considered, 
are input into Eq. (17), we obtain a residual value for 
each experiment, equal to the difference of the left- and 
right-hand sides of the equation. Thus, using Eq. 17, the 
25000 trials per value of 7 provide an ensemble of resid- 
ual values that are used to create a likelihood for 7 to be 
at that value, given the experimental inputs. The like- 
lihood, as a function of 7, can be obtained from x^(.l), 
where x^ = (m/''')^j M is the mean of the above ensem- 
ble of residual values, and a is the usual square root of 
the variance. The value of x^il) is then considered to 
represent a likelihood which is equal to that of a value x 
standard devations of a Gaussian distribution from the 
most likely value(s) of 7. We define the "exclusion level," 
as a function of the value of 7, as follows: the value of 7 
is excluded from a range at a given C.L. if the exclusion 
level in that range of 7 values is greater than the given 
C.L. 

We now turn to the VP decays. The method using 
VP decays shares the advantage with PP decays that no 
assumptions on strong phases are required. The disad- 
vantage is that, as we will see, the constraints from the 
VP modes are weak. 

We combine the information given above on the B" 
^ j)*±jjT^ B+ -^ D*+D^, and B+ -^ D+D*° branch- 
ing fractions and Acp information with measurements 
of the B" ^ D*-D+, B° -^ D- D*+ , B+ -> D*+W, 
and i3+ -^ DfD*'^ branching fractions [14], measure- 
ments of the B^ -^ D* D^ time-dependent CP asymme- 
tries [24, 37], and the world-average values of sin 2/3 [35] 
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and sin 0c [14]. Similar to the MC 7 determination for 
the VV and PP modes, we generate random values of 
each of the experimental inputs according to Gaussian 
distributions, with means and sigmas according to the 
measured central value and total errors on each experi- 
mental quantity. We again obtain a confidence level dis- 
tribution as a function of 7. 

Finally, we can combine information from the VV, PP, 
and VP modes. The resulting measured exclusion level 
as a function of 7 from each of the three sets of modes, 
as well as from their combination, is shown in Fig. 6. 
From the combined fit, we see that 7 is favored to lie in 
the range [0.07— 2.77] radians (with a +0 or -^-^T radians 
ambiguity) at 68% confidence level. This corresponds to 
[4.1°- 158.6°] (+0° or 180°). 

These constraints arc generally weaker than those 
found in Rcf. [8] due to the fact that the measured CP 
asymmetry in B^ -^ D*^D*^ has moved closer to the 
world-average sin2/3, with the newer i?° -^ D*'^D*^ 
measurements in Ref. [38]. The closer this CP asym- 
metry is to sin2/3, the weaker the resulting constraints 
are on 7, due to the fact that the closeness of the CP 
asymmetry to sin2/3 favors the dominance of the tree am- 
plitude, rather than the penguin amplitude whose phase 
provides the sensitivity to 7. Although the constraints 
are not strong, they contribute to the growing amount of 
information available on 7 from various sources. 



X. CONCLUSIONS 

In summary, we have measured branching fractions, 
upper limits, and charge asymmetries for all B me- 
son decays to D^*^D^*\ The results are shown in Ta- 
ble VII. This includes observation of the decay modes 
B" -^ D+D- and_B+ -^ D*+D*°, evidence for the decay 
modes B+ ^ D+D*^ and B+ -^ D+D" at 3.8ct and 4.9cr 
levels respectively, constraints on CP-violating charge 
asymmetries in the four decay modes B+ -^ £)(*)+ £)(*P ^ 
measurements of (and upper limits for) the decay modes 
B° -^ D*°D" and B° -^ D°D°, and improved branch- 



ing fractions, upper limits, and charge asymmetries in all 
other B -^ £)(*)£)(*) niodes. The results are consistent 
with theoretical expectation and (when available) previ- 
ous measurements. When we combine information from 
time-dependent CP asymmetries in B'^ -^ ]j{*)+]j{*)- 
decays [38, 39] and world-averaged branching fractions 
of B decays to Dl* D'-*' modes [14] using the technique 
proposed in Ref. [7] and implemented in Ref. [8] , we find 
the CKM phase 7 is favored to lie in the range [0.07—2.77] 
radians (with a +0 or +tt radians ambiguity) at 68% con- 
fidence level. 
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